Measurement of correlations between low-frequency vibrational modes and particle 
rearrangements in quasi-two-dimensional colloidal glasses 
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We investigate correlations between low-frequency vibrational modes and rearrangements in two- 
dimensional colloidal glasses composed of thermosensitive microgel particles which readily permit 
variation of sample packing fraction. At each packing fraction, the particle displacement covariance 
matrix is measured and used to extract the vibrational spectrum of the "shadow" colloidal glass 
(i.e., the particle network with the same geometry and interactions as the sample colloid but absent 
damping). Rearrangements are induced by successive, small reductions in packing fraction. The 
experimental results suggest that low-frequency quasi-localized phonon modes in colloidal glasses, 
i.e., modes that present low energy barriers for system rearrangements, are spatially correlated with 
rearrangements in this thermal system. 

PACS numbers: 82.70.Dd, 63.50.Lm 



In crystalline solids, plastic deformation is well under- 
stood to be associated with structural defects, partic- 
ularly dislocations lj. In amorphous solids, however, 
geometric structural parameters are not obviously con- 
nected to regions that later rearrange [2], hampering a 
priori identification of fragile regions. Recent numeri- 
cal simulations suggest that the spatial distribution of 
low frequency phonon modes may be correlated with 
irreversible rearrangements in glasses [3J 0], and that 
the quasi-localized low frequency vibrational modes of- 
ten observed in glasses [5] play a role in glass mechani- 
cal response [BHE]. In a related vein, experiments in vi- 
brated granular packings find that when cracks begin to 
appear, particles are likely to move in the direction of 
the polarization vectors of the lowest frequency modes 
4 . Taken together, these results suggest connections 
between rearrangement events and low-frequency vibra- 
tional modes in some glassy systems. They raise the 
question of whether such connections also exist in atomic 
and molecular glasses.. 

In this paper, we experimentally investigate corre- 
lations between vibrational modes and particle rear- 
rangements in quasi- two-dimensional (quasi-2D) colloidal 
glasses. The work is made possible by recently developed 
covariance matrix techniques that have led to the mea- 
surement of vibrational modes in thermal colloids from 
microscopy experiments [9llll|. The colloidal glasses 
are composed of thermosensitive microgel particles which 
readily permit in situ variation of sample packing frac- 
tion. At each packing fraction, the particle displacement 
covariance matrix of the sample is measured and used 
to extract the vibration spectrum of the "shadow" col- 
loidal glass, i.e., the particle network with the same ge- 
ometry and interactions as the sample colloid but absent 



damping. Particle rearrangements are observed to oc- 
cur when the packing fraction, or compressive stress, is 
decreased slightly. The experiments thus permit com- 
parison of packing-fraction-induced rearranging particle 
clusters with the spatial distribution of glassy phonon 
modes. Note that if our highly damped system were 
athermal, then rearrangements would only occur when 
an energy barrier-and, correspondingly, the frequency 
of a related vibrational mode-vanishes. In that case, 
the rearrangements would be deterministic. Our ther- 
mal system, however, can surmount low energy barriers 
associated with modes at nonzero frequency, leading to 
a stochastic response. Nonetheless, we find spatial corre- 
lations between regions of high displacement in low fre- 
quency quasi-localized phonon modes ("soft spots" [T2] ) 
and rearranging clusters of particles in glasses. Thus, 
our experiments illustrate that soft spots are correlated 
with rearrangements even when thermal fluctuations are 
present. Such resilience is critical if soft spots are to 
be viewed as the "structural defects" in disordered sys- 
tems, analogous to dislocations in crystals, that control 
mechanical response to applied load. 

The samples are disordered quasi-2D glasses com- 
posed of poly(N-isopropyl acrylamide) microgel colloidal 
spheres (i.e. NIPA particles) suspended in water. The 
NIPA particle diameter increases when temperature is 
reduced [13j [14] . Therefore, the sample packing fraction 
can be tuned over a relatively wide range, with suspen- 
sion evolving from colloidal fluid to jammed solid, via 
temperature changes of just a few degrees. The inter- 
particle potentials are measured to be short-range repul- 
sive with a soft tail |15) . enabling particles to be packed 
at high concentration. A binary mixture of microgel par- 
ticles, with diameters of 1 fan and 1.4 /xm and with a 
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FIG. 1: (color online) Particle displacement distribution when 
packing fraction changes from 0.885 to 0.879. Inset a: real 
space distribution of particle displacements; brightness of a 
particles is proportional to the displacement magnitude. Inset 
b: real space distribution of change of of the same region 
as Inset a, color change from black to red represents a change 
of ^6 from to 0.89, respectively. 

large/small particle number ratio of 0.7, was loaded be- 
tween two glass cover slips to create a colloidal mono- 
layer. The sample was hermetically sealed using optical 
glue (Norland 63) and annealed for 2 hours at 28 °C. 
Images of the samples were acquired at 30 frames/s us- 
ing bright field video microscopy at temperatures ranging 
from 24.9 to 26.3 °C. The sample temperature controls 
the proximity of the sample packing fraction to the jam- 
ming point, which is determined by the maximum of the 
height of the first peak in the pair correlation function 
to be <pj = 0.84 [15] . When the temperature is increased 
slightly (~0.2 °C), the particles shrink in size in less than 
0.1 s[16j. and they rearrange to find new equilibrium posi- 
tions [T7] as a result of unbalanced stresses in the system. 
A small decrease in particle size thus induces a uniform 
and unbiased perturbation of the system, similar to ther- 
mal annealing in bulk atomic and molecular materials. 
Sample temperature was controlled by thermal coupling 
to a resistively heated microscope objective (BiOptechs) 
and the sample was permitted to equilibrate for 15 min 
at each temperature before data acquisition. The ob- 
servation field-of-view was far from any boundaries or 
non-uniform regions within the sample, and trajectories 
of particles in the field-of-view were extracted using stan- 
dard particle tracking techniques (18) . Cage rearrange- 
ments were not observed during the 1000 second runs at 
constant packing fraction. Thus, the system remained 
in the same basin of the energy landscape. Brief rear- 
rangement events were observed as packing fraction was 
reduced, until the system settled down into another well- 
defined average configuration. 

The particle rearrangements by compression are spa- 



tially heterogeneous: some particles move, often in clus- 
ters, but most particles do not move measurable dis- 
tances. We resolve particle positions with 7 nm spa- 
tial resolution. The distribution of particle displacements 
when the packing fraction <fi was changed from tf> = 0.885 
to <j> = 0.879 is plotted in Fig. 1. Notice that most par- 
ticles experience a displacement of less than 40 nm, i.e., 
less than 1/30 of a particle diameter. In our analysis, 
these particles are considered to have remained in their 
equilibrium position during the temperature change; this 
non-zero displacement is due to the drifting of equilib- 
rium positions as a result of elastic relaxation, and the 
diffusion of particles within cages over long periods of 
time (~2000 seconds). Such drifting is typically less than 
40 nm, and it does not amount to a structural rearrange- 
ment. These drifts, therefore, do not influence the struc- 
tural or dynamical relation between particles, since the 
particles remain in the same cages of their neighbors. 
The distribution has a long tail due to particles which 
have experienced significant displacements. The latter 
particles tend to be spatially clustered as shown in Fig. 1 
(Inset a). We also measure the change in the local bond 
orientational order parameter, i.e., for the jth particle, 
^6 =w~ Y^k nn e 6l9jk , summed over all nearest neighbor 
bonds; 9jk is the angle between the bond joining parti- 
cles j, k and the x-axis, and N nn is the number of near- 
est neighbors. In Inset b of Fig.l, the distribution of the 
absolute value of each particle's change in bond orienta- 
tional order, is plotted. This peak in the change of bond 
orientational order is spatially correlated with clusters of 
particles with large displacement. Thus, both particle 
displacements and changes in bond orientational order 
demonstrate the heterogeneous nature of rearrangements 
in the packing. 

The measured displacement covariance matrix [9]- 
HT| fT9] was employed to extract vibrational modes of the 
"shadow" colloidal system, which shares the same geo- 
metric configuration and interactions of the experimental 
colloidal system, but is undamped. A displacement co- 
variance matrix is constructed by calculating the equal- 
time covariance of displacements between each particle 
pair for each direction. The relative amplitude and di- 
rection of motion due to a mode u> is thus described in 
terms of the polarization vector ef for each particle i. 

Fig. 2 exhibits colored contour plots of polarization 
magnitudes for each particle summed over the 30 lowest 
frequency modes, i.e., ^I]^°(e" 7 ) 2 , where is the 7 
coordinate of the polarization vector e*i of particle i for 
a normalized eigenmode with frequency lj, as in Ref. 3. 
These vibrational properties correspond to the packing 
configuration before the packing fraction change. In ad- 
dition, the circles plotted in the figure correspond to 
particles that experienced significant displacements (i.e., 
displacements larger than 40 nm) following the packing 
fraction change. Use of the change in gives a sim- 
ilar selection of clusters. Since we are most interested 
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FIG. 2: (color online) Rearranging clusters and low-frequency modes at three different packing fractions studied, ranging from 
4> = 0.885 to (j> = 0.850 from a to c. Color contour plots indicate polarization magnitudes for each particle, summed over the 
lowest 30 modes. Circles indicate particles that moved more than 40 nm at each decompression step. 



in large rearrangement events, we group these displaced 
particles into clusters based on nearest neighbor pair- 
ings, with clusters containing fewer than five particles 
ignored. A good correlation between the spatial distri- 
bution of the low frequency modes and the rearranging 
clusters is apparent over several steps of change in pack- 
ing fraction [17]. (Note, these results are not sensitive 
to the number of modes chosen; similar results are ob- 
tained using the lowest 10 or 50 modes.) We see that the 
number of rearranging clusters, as well as the total num- 
ber of rearranging particles, increases as </> is decreased 
towards the jamming transition as the system gradually 
loses its mechanical stability. Note that participation in 
low frequency modes correlates with regions of substan- 
tial particle rearrangement, as found in Ref. 3. How- 
ever, rearrangements studied in Ref. 3 only involved par- 
ticles that irreversibly changed four nearest neighbors (or 
more), and involved samples below the jamming point. 
In our experiments, no particle experienced more than 
four nearest neighbor changes. 

The correlation between the vibrational modes and re- 
arrangements can be quantified by the normalized inner 
product between the eigenvectors, which forms a com- 
plete basis for 2N-component vectors [4], and the dis- 
placement vector R c i that contains the displacement com- 
ponents of all particles in the field of view. Here N is the 
number of particles in the field of view. A projection 

coefficient is defined as a 2 (uj) — , where is 

\R 2 J 

the eigenvector of mode u. a 2 (oj) for all packing frac- 
tion steps is plotted in Fig. 3. Clearly, the contributions 
of the lowest frequency modes (particularly the lowest 
N m R3 10 modes) to the rearrangement vectors are much 
larger than those from higher frequency modes, and they 
are more than two orders of magnitude larger than the 
correlation noise of 1/2N. Furthermore, the accumulated 
projection coefficients of the lowest frequency modes are 
significant [17]. Thus, rearrangements are not only likely 
to occur in peak regions of low frequency modes, but 
the displacements of individual particles also tend to fol- 
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FIG. 3: (color online) Projection coefficient a 2 versus mode 
number. Different symbols represent a 2 from different pack- 
ing fraction steps. The eigenvectors are derived at packing 
fractions as indicated by the symbols, and the rearrangement 
vector is determined from the displacements of particles when 
the system is decompressed from the packing fraction indi- 
cated to the next highest value. The data are binned with a 
window of 5. 

low the directions of the polarization vectors of the these 
modes, consistent with numerical findings for sheared 
systems [7j. This experimental finding strongly supports 
the notion that the quasi-localized low frequency modes 
present the lowest energy barriers to collective motion 
0- 

Thus far, our experiments suggest that rearrangements 
in disordered colloidal packings tend to occur in regions 
where low frequency phonon modes concentrate. How- 
ever, not all such regions experience rearrangements, as 
shown in the high packing fraction results in Fig. 2a and 
b. Recent numerical work on athermal packings suggests 
that the mechanical response may be controlled by a pop- 
ulation of soft spots, consistent with phenomenological 
theories of shear transformation zones [20] , where at least 
one soft spot is correlated to a given rearrangement |12j . 
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FIG. 4: (color online) Comparison of soft spots and rearranging clusters at high packing fractions. Blue solid circles indicate 
particles belonging to rearranging clusters, defined as particles with displacements exceeding 40nm; red empty circles are 
particles belonging to soft spots for N m = 10, N p = 130. 



We now show that soft spots also exist in our colloidal 
system at finite temperatures. 

Simulations in Ref [12] show that regions prone to re- 
arrangement or "soft spots" are characterized by an im- 
portant energy scale and length scale. The energy scale 
is characterized by the vibrational modes that are most 
relevant to particle rearrangements, and the length scale 
is associated with the number of particles participating in 
a rearranging event. Following Ref. [12], Fig. 3 suggests 
that the lowest frequency vibrational modes, particularly 
the lowest N m — 10 of them, are most correlated with 
particle rearrangements; thus, the frequency of the 10th 
lowest frequency mode provides an estimate of the energy 
scale. Similarly, the experiments show that the largest re- 
arrangement observed contains N p ss 130 particles; this 
determines the characteristic length scale. These choices 
for N m and N p are not arbitrary but have physical moti- 
vation; we also note that the results are quite insensitive 
to the particular values of N m and N p , so that similar 
results are obtained for N m = 50 and N p = 65 [17] • We 
limit our analysis to the high packing fraction regime in 
order to best match the conditions in simulations wherein 
only one rearrangement event occurs at a time. 

The resulting soft spots are shown in comparison with 
the rearranging clusters in Fig. 4. At high packing frac- 
tions, rearrangements reliably occur at soft spots, sug- 
gesting that soft spots are pre-existing structural defects 
wherein the material preferentially fails. Note also that 
the average soft spot size is roughly 20 particles and the 
density of the soft spots is about 10% of the total num- 
ber of particles in the field of view at all packing fractions 
studied. 

Fig. 4 shows that while the soft spot analysis identifies 
a population of regions wherein the particles are likely to 
rearrange, it does not predict which of those regions will 
rearrange. In a thermal system such as ours, the precise 
location of the rearrangement is chosen stochastically, 
depending on the energy barriers encountered. This sce- 
nario differs fundamentally from the situation for systems 
sheared quasistatically at zero temperature, in which soft 
spots were previously identified |12j . In addition, we note 



that the simulation focused only on the first initial rear- 
rangement event with no regard to possible avalanches 
or separate excitations at slightly later times. Experi- 
ments, on the other hand, are susceptible to avalanches 
and necessarily consider all rearrangement events, includ- 
ing secondary motions, that occur during the change of 
packing fraction. Nonetheless, Fig. 4 suggests that suffi- 
ciently far above the jamming transition, rearrangements 
tend to occur at soft spots even in our thermal system, 
consistent with the idea that soft spots are robust, intrin- 
sic defects in disordered solids that control where local 
failure occurs under mechanical load. The experiments 
thus raise the possibility of using vibrational properties 
to identify regions susceptible to failure in packings of 
constituents ranging from nanoparticles to macroscopic 
grains, perhaps even including systems subject to natural 
events such as earthquakes and avalanches [2D H2] • 
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